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Abstract. We determine ages of 71 old Open Clusters by a two-step method: we use main-squence fitting to 10 
selected clusters, in order to obtain their distances, and derive their ages from comparison with our own isochrones 
used before for Globular Clusters. We then calibrate the morphological age indicator 5(V), which can be obtained 
for all remaining clusters, in terms of age and metallicity. Particular care is taken to ensure consistency in the 
whole procedure. The resulting Open Cluster ages connect well to our previous Globular Cluster results. From 
^ ' the Open Cluster sample, as well as from the combined sample, questions regarding the formation process of 

Galactic components are addressed. The age of the oldest open clusters (NGC 6791 and Be 17) is of the order of 
■ 10 Gyr. We determine a delay by 2. Oil. 5 Gyr between the start of the halo and thin disk formation, whereas thin 

and thick disk started to form approximately at the same time. We do not find any significant age-metallicity 
relationship for the open cluster sample. The cumulative age distribution of the whole open cluster sample shows 
a moderately significant (~ 2a level) departure from the predictions for an exponentially declining dissolution rate 
with timescale of 2.5 Gyr. The cumulative age distribution does not show any trend with galactocentric distance, 
but the clusters with larger height to the Galactic plane have an excess of objects between 2-4 and 6 Gyr with 
' respect to their counterpart closer to the plane of the Galaxy. 

q . Key words. Galaxy: disk - evolution - open clusters and associations: general - stellar content 
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' 1. Introduction absolute magnitudes), which have to be fitted by the ap- 

propriate theoretical model - stems from the fact that star 
The theory of the formation of galaxies is without any clusters are made of coeval objects, largely with the same 
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doubt one of the outstanding problems of astrophysics. mitial chemica i composition and located at the same dis- 

. Although in the past decades considerable progress has tanC6j so that it ig possible to use morphological parame- 

been made, we do not have yet a complete and defini- ters de duced from theoretical isochrones in order to derive 

tive picture of how galaxies form. As discussed by, e.g., their age In this way one can bypass the thorny prob i em 

Freeman & Bland-Hawthorn (2002), a detailed study of of determining a reliable empirical and theoretical temper- 

the formation of the Galaxy lies at the core of under- ature scalej and of acquiring high resolution spectroscopy 

standing the complex processes leading to the formation f Qr j ar g e sam pl e s of stars 
of external galaxy systems. A way to shed some light on 

this problem is to study the timescale for the formation of In a series of papers published in the last 6 years 

the different Galactic populations, e.g., halo, thick disk, (Salaris et al. 1997; Salaris & Weiss 1997, 1998; Salaris 

thin disk and bulge, by means of stellar age dating. The & Weiss 2002, hereinafter SW02), we have addressed the 

most reliable stellar ages are obtained for the star clus- problem of the timescale for the formation of the halo 

ters belonging to the various populations, i.e., the globu- and thick disk by homogeneously determining the age of 

lar clusters (GCs) in the halo, thick disk and bulge, and a large sample of Galactic GCs. The latest SW02 study 

the open clusters (OCs) in the thin disk. The advantage (including 55 GCs) concluded that metal poor clusters 

of dating star clusters over individual stars - whose age (up to [Fe/H] between -1.6 and -1.2, depending on the 

determination relies entirely on the knowledge of individ- adopted [Fe/H] scale) are coeval within -1 Gyr, with an 

ual metallicities, effective temperatures and gravities (or age of the order of 12-13 Gyr, whereas the more metal 

rich ones show an age spread, are on average younger and 

Send offprint requests to: M. Salaris display a weak age-metallicity relationship (age decreasing 
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with increasing [Fe/H]). This result is in agreement with 
other independent analyses, such as that by Rosenberg et 
al. (1999). When searching for relationships between age 
and position within the halo, it was found that the age 
spread starts from galactocentric distances (i? gc ) between 
8 and 13 kpc outwards, the precise value depending again 
on the adopted [Fe/H] scale. 

It is now important to address the question of when the 
thin disk started to build up, relative to the thick disk and 
halo. This can be accomplished by studying the age distri- 
bution of the oldest OCs. In general, OCs are expected to 
be disrupted easily by encounters with massive clouds in 
the disk (Spitzer 1958); however, the most massive OCs 
or those with orbits that keep them far away from the 
Galactic plane for most of their lifetimes are expected to 
survive for longer periods of time. These old objects are 
therefore test particles - in analogy to the GCs - prob- 
ing the earliest stages of the formation of the disk. It is 
essential to determine their ages using stellar models and 
methods which place them on the same scale as GC ages. 
An analysis of this kind, based on homogeneous age dat- 
ing of all the known old OCs and a large sample of GCs, 
employing the latest generation of stellar models is still 
lacking (see, e.g., Liu & Chaboyer 2000 for a study of this 
kind, but considering only a very small number of OCs 
and GCs), and this paper is intended to fill this gap. 

Here we will reanalyze the old OC sample reviewed 
by Friel (1995, hereafter F95), and based on the seminal 
papers by Phelps, Janes & Montgomery (1994) and Janes 
& Phelps (1994, hereafter JP94), to which we have added 
two additional clusters (ESO 093-SC08 and vdBH 176) 
recently studied by Phelps & Schick (2003). This should 
contain approximately all presently known old OCs. Our 
aim is to determine their age on a scale consistent with 
the GC ages determined by SW02, to study the existence 
of possible relationships between age, position within the 
disk and [Fe/H], and to compare their ages with the GC 
population. In Sect. 2 we describe the cluster sample and 
the techniques used to determine their age. The result- 
ing age distribution is analyzed in Sect. 3, while Sect. 4 
deals with the comparison with the GC ages by SW02. A 
summary and conclusions follow in Sect. 5. 

2. Cluster sample and age determination method 

We consider a total of 71 clusters - 69 from Friel (1995) 
and 2 from Phelps & Schick (2003) - whose morphologi- 
cal age parameter S(V) is equal to or larger than the value 
for Praesepe, i.e. S(V)=0.3 (see below for the definition of 
S(V)). The definition of old OCs by JP94 is slightly differ- 
ent: they considered as 'old' all OCs where 5(V) is larger 
than zero. We have not included objects with 8(V) <0.3, 
because we did not have clusters with 5(V) between and 
0.3 that could be used to calibrate adequately a relation- 
ship between this parameter and cluster age (see below 
for details about the calibration). 

Friel's (1995) sample is mainly the same as the one 
studied by JP94, with only a few additions. JP94 have 
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Fig. 1. Comparison between the distances given by JP94 
and the MS fitting distances we obtain for a subsample of 
10 clusters (see text for details). The solid line represent 
the 1:1 relationship between the two sets of distances. 



discussed in detail the completeness of their sample of old 
OCs, and concluded that most probably the number of 
undetected old clusters is small and should not have a 
major effect on the overall age distribution, even though 
the properties of the age distribution perpendicular to the 
Galactic plane may be affected by still undetected old 
OCs, which should be preferentially located very close to 
the plane of the Galaxy. 

The [Fe/H] values for our sample are taken when- 
ever possible (38 clusters) from Gratton (2000, hereafter 
GOO), who transformed various metallicity scales based on 
low resolution spectroscopy onto an homogeneous scale 
tied to high resolution [Fe/H] determinations. One ex- 
ception is the cluster Praesepe, for which we have em- 
ployed the Hyades metallicity (GOO reports a lower value), 
based on the discussion and references in Percival et 
al. (2002, hereafter PSK02). In case of clusters not listed 
by GOO we have either used the value provided by the 
WEBDA OC database | |http://obswww.unige.ch/webda/[ 
see Mcrmilliod 1992) when available, to which we at- 
tached an error of 0.15 dex (11 clusters), or we as- 
sumed [Fe/H] =0.0 with an error of 0.20 dex (22 clusters). 
The cluster galactocentric distances and heights to the 
Galactic plane are taken from Friel (1995) and Phelps & 
Schick (2003). For many of the clusters in our analysis the 
existing photometry and/or uncertainties in the cluster 
parameters do not allow to perform a more accurate and 
homogeneous distance determination. Therefore we used 
the results presented in the mentioned papers, where more 
details about this issue can be found. 
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As a test we compared in Fig. ^ the distances we ob- 
tained for a subsamplc of 10 clusters (see next subsection 
for details) from the Main Sequence (hereinafter MS) fit- 
ting technique, with the JP94 results. We did not find any 
statistically significant trend of the difference between the 
two sets of distances with respect to our MS fitting de- 
terminations. The mean value of the difference is equal to 
only 9 pc, with a dispersion of 160 pc around the mean. 

The data for the complete cluster sample are summa- 
rized in Table 1; the flag attached to each cluster pro- 
vides the source for the metallicity; a value equal to or 
1 means [Fe/H] from GOO, whereas a value of 2 means 
that the source is the WEBDA database or that there is 
no available [Fe/H] determination. A flag equal to de- 
notes the subsample of clusters that are used for our age 
calibration, as explained in the next subsections. 

We have used as age indicator the morphological pa- 
rameter 5(V) defined by JP94, which is similar to the 
A(V) parameter used in GC dating (e.g. SW02), cali- 
brated in terms of absolute age and [Fe/H] following the 
same kind of approach as in JP94, and in Carraro & 
Chiosi (1994a) for their analysis of a sample of 36 old 
OCs. 

2.1. The morphological age index 5(V) 

The use of morphological indices that quantify differences 
in the Colour-Magnitude Diagram (CMD) of clusters in 
terms of age differences is a well established technique 
(see, e.g., Anthony- Twarog & Twarog 1985 and JP94 for 
OCs; Rosenberg et al. 1999 for GCs); it allows to establish 
a relative age ranking among a given cluster sample, by- 
passing the well known difficulties with isochrone fitting 
methods (see, e.g., VandcnBcrg et al. 1990; Sarajedini & 
Demarque 1990; Salaris & Weiss 1997, SW02). 

Phelps et al. (1994) and JP94 have defined two mor- 
phological parameters, called <5(1^) 1 and 51 which they 
applied to their sample of old OCs. 5(V) is defined as 
the magnitude difference between the cluster turn-off re- 
gion and the He-burning clump stars. More precisely, the 
reference point in the turn-off region is taken as the inflec- 
tion point between the turn-off and the base of the giant 
branch. This point is well defined and unaffected by the 
presence of a binary sequence and/or field stars, accord- 
ing to Phelps et al. (1994). 51 is the difference in colour 
index between the bluest point on the MS at the luminos- 
ity of the turn-off and the colour of the giant branch one 
magnitude brighter than the turn off luminosity. A sim- 
ple linear relationship between 5(V) and 51 was found by 
JP94 when analysing the clusters where both indices could 
be measured, and this was applied by the same authors to 
obtain an estimate for 5(V) for those clusters where clump 
stars were not identified. This estimate we will call 5(V) 
in the following. JP94 choose 8(V) as the primary age 
indicator, so that for clusters without visible clump the 
estimated value described above was used; for the other 

1 In the quoted paper this quantity is denoted as 5(V) 



clusters the final 5(V) given by JP94 is the average be- 
tween the observed 5(V) and the one computed from the 
51 — 5(V) conversion described above. 

The 6(V) data for our cluster sample are reported in 
Table 1; the associated errors are derived from the quality 
grade assigned by JP94. Following JP94 we considered 
errors by, respectively, 0.05 mag for clusters graded "a" , 
0.15 mag for clusters graded "b" and 0.25 mag for clusters 
graded "c". In case of the two clusters from Phelps & 
Schick (2003) we have assumed an error by 0.25 mag. 

2.2. Calibration of 5(V) 

The cluster 5(V) values given in Table 1 can be translated 
into absolute ages by determining a relationship 5(V)-t- 
[Fe/H], based on a subsample of clusters with high qual- 
ity CMDs, spanning the entire [Fe/H] and 5(V) range of 
the full cluster sample, and for which the age can be de- 
termined with confidence. JP94 determined a relationship 
between 5{V) and age - neglecting the effect of metallicity 
- based on a sample of OCs and GCs with age determi- 
nations obtained by various authors and with a variety of 
methods and stellar models. In case of multiple age deter- 
minations for the same objects JP94 averaged the results 
from the various authors. As clearly stated by JP94, due 
to the heterogeneity of the calibration material, their cali- 
bration was mainly aimed at producing the ranking of the 
clusters in terms of relative ages. 

Here we wish to obtain a new highly homogeneous and 
reliable calibration in terms of absolute ages, based on a 
consistent set of updated stellar models. We have consid- 
ered a subsample of 10 OCs (clusters with the flag value 
equal to in Table 1) plus 1 GC (47 Tuc), whose ages 
have been determined by fitting the CMD turn-off lumi- 
nosity with theoretical isochrones, after determining their 
distance from an empirical MS fitting technique which 
employs large samples of field MS stars with accurate 
Hipparcos parallaxes. In this way the ages we obtain are 
firmly tied to the Hipparcos distance scale. 

The stellar models used to determine the OC ages have 
been computed with exactly the same updated physics 
employed for calculating the GC isochrones by Salaris 
& Weiss (1998); these GC isochrones have been used to 
determine the ages of the large GC sample analyzed by 
SW02, and we refer the reader to Salaris & Weiss (1998) 
for details about the model input physics. The turn-off 
stars in the younger clusters in our sample do have con- 
vective cores, and therefore we have included in our models 
overshooting beyond the formal boundary of the convec- 
tive core (i.e., instantaneous mixing and radiative tem- 
perature gradient in the overshooting region beyond the 
boundary of the convective core fixed by the Schwarzschild 
criterion), with an extension of 0.2 pressure scale heights 
for masses above 1.4M Q , and linearly decreasing to zero 
from M=1.4M to M=1.0M Q . This prescription is in 
broad agreement with the results obtained by Ribas et 
al. (2000) from the comparison of stellar models with 
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eclipsing binary systems, at least in the mass interval 
spanned by the stars evolving in the turn-off region of 
our OC sample. We have computed stellar models and 
isochrones for the appropriate metallicity (scaled solar 
Grevesse & Noels 1993 heavy element distribution) of each 
of our calibrating clusters, using an He-mass fraction Y 
that follows the relationship Y = 0.248 + 1.44 Z. The pri- 
mordial He is derived from the recent results of analyses 
of the CMB power spectrum (see, e.g., the discussion in 
Cassisi et al. 2003), whereas the slope A(Y)/A(Z) = 1.44 
arises from the constraint imposed by the initial He- 
abundance of the standard solar model. In case of 47 Tuc 
and the other GC ages discussed in Sect. 4 we have com- 
puted selected a-enhanced isochrones (the same metal dis- 
tribution as in Salaris & Weiss 1998) at various metal- 
licities and with the same Y—Z relationship as for the 
OCs, in order to determine the age of 47 Tuc and revise 
the ages of the GCs in SW02 which were computed using 
Y = 0.230 + 3.0 Z. 

The MS fitting distances to the calibrating OCs are 
from PSK02 and Pcrcival & Salaris (2003, hereafter PS03), 
with the exception of NGC 6791 (see below), whilst the 
distance to 47 Tuc is from Percival et al. (2001, hereafter 
P01); they are all based on two large samples of unevolvcd 
field MS stars with accurate Hipparcos parallaxes and in- 
dividual metallicity determinations. The more metal rich 
sample of 54 dwarfs with [Fe/H] between ~ —0.4 and 
~ +0.3 (field dwarfs) has been used to derive the dis- 
tances to the calibrating OCs, while the distance to 47 Tuc 
has been obtained using a sample of 43 more metal-poor 
dwarfs with [Fe/H] between ~ —1.0 and ~ —0.3 (sub- 
dwarfs). Table 2 contains the results from the MS fit- 
ting distance determinations; reddenings are from Twarog 
et al. (1997) with an associated error of ±0.02 mag as 
adopted by Sarajedini (1999). It is important to notice 
that, as discussed in PSK02, our MS fitting distances to 
the Hyades and Praesepe agree well with the Hipparcos 
parallax measurements. 

A detailed description of the distance determination 
method is given in the three papers mentioned above, to- 
gether with the sources for the adopted MS CMDs. We 
just recall here that the method is based on constructing 
an empirical template MS from the field stars, by apply- 
ing colour shifts to the individual objects to account for 
the differences in metallicity between the field stars and 
the cluster. The template is then shifted in magnitude 
to match the dereddened and extinction-corrected clus- 
ter MS, the extent of the shift being equal to the distance 
modulus (m—M)o- The colour shifts for the field star sam- 
ple have been derived empirically as discussed in PSK02 
and, strictly speaking, they are applicable only within the 
metallicity range spanned by the field dwarfs themselves, 
i.e., up to [Fe/H]~0.3. In case of the subdwarf sample 
used for 47 Tuc, the shifts have been obtained from the 
differential use of the isochrones by Salaris & Weiss (1998) 
as discussed in P01. Whenever possible we have derived 
MS fitting distances using both the (B — V) and (V — I) 
colours, always finding agreement (within the associated 
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Fig. 2. Best-fitting isochrone for the MS and turn-off re- 
gion of the cluster NGC 2420. Reddening, distance mod- 
ulus, age and metallicity employed in the fit are given. 

error bars, which, for given reddening and [Fe/H] values, 
are typically of a few hundredths of magnitude, and usu- 
ally the same for both colours) between the two values. 
In this case the final distance is the unweighted mean be- 
tween the results in (B — V) and in (V — I). The final 
error budget takes into account the errors in the cluster 
reddening and [Fc/H] quoted in Tables 1 and 2. 

P01 and PSK02 discussed at length the consistency be- 
tween the metallicity scale of the field stars and clusters for 
each of the two separate samples, a necessary prerequisite 
for the reliability of the MS fitting distances. What mat- 
ters most in our case is the consistency between the dis- 
tances obtained separately for 47 Tuc with the subdwarf 
sample, and for the other calibrating OCs with the field 
star sample. In fact, they have been determined from two 
different samples of field MS stars, with metallicity scales 
determined independently, and using different methods to 
determine the colour shifts. Consistency between the two 
sets of distances means that one should be able to use 
the subdwarf sample to derive the MS fitting distance 
modulus to a cluster like, e.g., NGC 2420, which is at 
the lower metallicity end of the calibrating OCs, and re- 
cover the value obtained employing the field dwarf sam- 
ple (11.94+0.07 mag). Since we get consistent distances 
in both (B - V) and (V - I) for the distance to 47 Tuc 
and to other calibrating OCs, we performed this test with 
the (V — I) colour only. If consistency is achieved for this 
colour, it is automatically ensured for (B — V), too. 

We performed this test using the entire subdwarf sam- 
ple and determining the appropriate colour shifts from 
the theoretical isochrones discussed before. We finally ob- 
tained a distance modulus that is within 0.02 mag of the 
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circles) given in Table 1. 



value obtained from the field dwarf sample, confirming the 
consistency between the distances obtained from the two 
separate samples of subdwarfs and field dwarfs. 

We also included NGC 6791 in the calibrating sample, 
in order to extend our calibration to very high metallici- 
ties, and have another extremely 'old' (e.g. JP94) calibrat- 
ing cluster in addition to 47 Tuc. We used E(B — V) = 
0.15±0.02 (Twaroget al. 1997), [Fe/H]=0.40±0.06 (GOO), 
and employed the same field dwarf sample and method as 
for the other OCs; the cluster MS CMD is from the new 
photometry by Stetson et al. (2003). We determined the 
distance modulus using both the (B — V) and (V — I) 
colours (they provide the same distance modulus within 
0.01 mag), obtaining (m — M)o=12.96±0.10, as reported 
in Table 2. 

As a note of caution we notice that the cluster 
NGC 6791 metallicity is slightly above the upper limit 
of the [Fe/H] range where our MS fitting method is appli- 
cable, so that we had to slightly extrapolate the empirical 
colour shifts applied to the field dwarfs. However, we have 
tested that the exclusion of NGC 6791 from the calibrat- 
ing sample does not alter substantially the calibration of 
the sought 5(V)-t-[Fe/H] relationship and therefore we re- 
tained NGC 6791 in our calibrating sample. 

After the MS fitting distances have been determined, 
cluster ages for the calibrating clusters have been obtained 
from isochrone fitting to the CMD turn off region, and are 
given in Table 1. Figure [3 shows an example (the clus- 
ter NGC 2420) of our age determination. The error bar 
for the age includes in quadrature the contributions due 
to the uncertainty in the cluster distance modulus and 
metallicity. 



We notice that our ages are very similar to the values 
obtained by PS03 using the same distance moduli but the 
isochrones by Girardi et al. (2000). We also, as a test, 
determined the ages of our calibrating clusters by means 
of the Lejeune & Schaerer (2001) isochrones, obtaining 
ages within less than 10% of the values given in Table 1. 

With ages, [Fc/H] (all on the GOO scale) and 6(V) val- 
ues of the calibrating clusters we determined the sought 
calibration of age as a function of 8(V) and metallicity. 
We found a simple relationship (shown in Fig. |3J) between 
the logarithm of the cluster age and both [Fe/H] and S(V), 
given by 

log(t) = 0.04 5{Vf + 0.34 5(V) + 0.07 [Fc/H] + 8.76 (1) 

with a lcr dispersion equal to 0.062 dex. 

The dependence of the logarithm of the age on S(V) is 
not very different from the calibration by JP94, who found 
log(i) to be proportional to 0.256 S(V) + 0.0662 5{V) 2 , 
without including a metallicity term. Interestingly, 
Carraro & Chiosi (1994a) calibrated log(i) in terms of a 
morphological parameter similar to S(V), and determined 
a dependence on metallicity equal to 0.08 [Fe/H], almost 
identical to our result for the JP94 S(V). 

It is interesting to notice that the [Fe/H] dependence of 
Eq. (1) is qualitatively the same as in theoretical models. 
In fact, by using, e.g., the Girardi et al. (2000) isochrones, 
we have computed the magnitude difference between He- 
burning clump and turn off for the age and [Fe/H] range 
of the studied clusters. We found that, for a given value of 
this magnitude difference, metal poorer isochrones provide 
lower ages, as predicted by Eq. (1). The reason is that in 
this [Fe/H] range, once the age is fixed, a decrease of the 
metallicity increases the clump luminosity more than the 
turn off one. Owing to the fact that the clump level is 
weakly dependent on age for our relevant age range, it 
is therefore clear that isochrones with lower metallicity 
have to be younger in order to show the same clump- 
turn off magnitude difference as more metal rich ones. The 
opposite happens in the regime of metal poor globular 
clusters, where the dependence of the horizontal branch 
magnitude on metallicity is weaker than the turn off one 
(at fixed age). 

3. The age of the old OCs 

In order to determine the ages of the remaining 61 OCs 
in our sample we have applied Eq. (1) to their S(V) and 
[Fe/H] values displayed in Table 1. 

The error bar on the individual determination has been 
computed by combining in quadrature the contribution to 
the uncertainty arising from the dispersion associated to 
Eq. (1), plus the contribution due to the [Fe/H] and 8(V) 
error propagation through Eq. (1). The uncertainty in the 
metallicities does not play a significant role in the error 
budget, due to the weak dependence of log(i) on [Fe/H]. 

In Fig.0]we plot the cluster ages against their [Fe/H]. 
Open squares denote the 11 calibrating clusters (10 OCs 
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Table 1. Cluster data. The columns display, respectively, cluster name, value of the 5(V) morphological parameter 
and its associated error, [Fe/H] and associated error, age in Gyr and associated error, galactocentric distance in kpc, 
height to the Galactic plane in pc, source of [Fe/H] value, previous age estimate on the JP94 scale (see text for details). 
The last 11 clusters are the calibrating clusters for our t-[Fe/H]-<5(V) relationship. 
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plus 47 Tuc) which cover the entire [Fe/H] and age range 
spanned by the full sample, with the only exception of 
Be 21, which has a metallicity [Fe/H] ~ —1.00, and is 
therefore outside the range of validity of Eq. (1); thus its 
age has to be treated with caution. 

Be 17 and NGC 6791 appear to be the oldest known 
OCs, and their ages are formally the same, within the error 
bars, as the age of the thick disk GC 47 Tuc (see Sect. 4 for 
a discussion about the comparison with GC ages) . FigureEl 
compares our ages with the results by JP94. Our values 
are, with few exceptions, systematically lower, with the 
oldest clusters - Be 17 and NGC 6791 - having ages of 
about 10 Gyr, whereas the oldest cluster in JP94 is Be 17, 
with an estimated age of 12.6 Gyr, about 2.5 Gyr higher 
than our results; our age for Be 17 (10.1 ± 2.8 Gyr) is in 
line with the recent analysis by Carraro et al. (1999) who 



determined a value of 9 ± 1 Gyr from isochrone fitting, 
although our estimate (based on the S(V) provided by 
JP94) has a much larger error bar. Our age for NGC 6791 
- 10.2 ± 1.2 Gyr - is larger than the result by Chaboyer et 
al. (1999) who found 8.0 ± 0.5 Gyr from isochrone fitting, 
and compatible, within the error bar, with the age of 8- 
9 Gyr determined by Carraro et al. (1999). 

The lower panel of Fig.Elshows a histogram of the clus- 
ter ages (without the GC 47 Tuc, that in the following will 
always be excluded from the analysis of the OC sample). 
As mentioned before, according to JP94 the shape of this 
distribution for the whole sample should not be altered 
by the still undetected old OCs. The corresponding ob- 
served cumulative function (i.e., number of clusters with 
ages larger than a given value U) is displayed in the upper 
panel of the same figure (open squares). 
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Fig. 4. Ages for the 71 old OCs plus 47 Tuc, as a function 
of the cluster [Fe/H]. Open squares denote the 11 calibrat- 
ing clusters. Individual clusters discussed in the text are 
labelled. 



These data contain in principle important information 
about the timescales of cluster destruction. In the sim- 
plest case of a uniform formation rate and exponentially 
declining dissolution rate, the open squares in the upper 
panel of Fig. [S] should follow a straight line whose slope 
is equal to the inverse of the dissolution timescale (see, 
e.g., Janes et al. 1988). In fact, the situation may be more 
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t jp94 ( G y r ) 

Fig. 5. Comparison between our derived ages and those 
given by JP94. The solid line denotes the 1:1 relationship. 

complex, for the points in the cumulative age distribution 
apparently do not follow a single slope; there is a change 
in the shape of the cumulative age distribution, visible 
between ~4 and ~6 Gyr, which appears as a 'bump' in 
the differential age distribution. This 'excess' of clusters 
has been already noticed by JP94, and it was located in 
the age interval between 5 and 7 Gyr on their age scale. 
We have investigated further this matter by evaluating the 
error bars associated to the points in our cumulative age 
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Fig. 6. Cumulative age distribution (upper panel) and dif- 
ferential age distribution (lower panel) for the 71 OCs in 
our sample. In the upper panel the age distribution from 
the Monte Carlo simulations discussed in the text (filled 
circles) is shown together with the actual cumulative dis- 
tribution (open squares). The line with slope correspond- 
ing to a dissolution timescale of 2.5 Gyr is also displayed. 



distribution. For this purpose, we performed an extensive 
Monte Carlo simulation by considering the individual OC 
ages given in Table 1, together with their associated er- 
rors. We then determined 10000 synthetic samples of ages 
for our 71 clusters. In each sample the individual clus- 
ter ages were randomly assigned according to a Gaussian 
distribution centred around the values given in Table 1, 
with a la dispersion equal to the estimated individual er- 
rors. We then determined the cumulative age distribution 
for each of the 10000 synthetic samples - using the same 
age bins as in Fig. 0- and determined the mean number 
counts and associated la dispersion in each age bin (an 
analogous result is obtained if we use the modal value of 
the number counts in each bin). 

The upper panel of Fig. displays the resulting syn- 
thetic cumulative age distribution as filled circles. It is 
worth noticing that in general the actual distribution lies 
comfortably within the la error bars associated to the 
synthetic one, as expected, since it has to correspond to 
one realization of our ensemble of synthetic age distribu- 
tions. We have also displayed the best fit single slope that 
matches the data: it corresponds to a dissolution timescale 
of 2.5 Gyr. Once the effect of the age errors is taken into 
account, we used a % 2 test to determine the significance 
of the excess of clusters in the age range between 4 and 
6 Gyr, which results to be at ~ 2a level for the whole 
sample of objects. 



Table 2. Distance moduli and reddenings of the calibrat- 
ing clusters (see text for details). 
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Fig. 7. Comparison of the cumulative age distributions 
for clusters located in two different ranges of height above 
the Galactic plane, \z\ (upper panel), or galactocentric 
distance, R gc (lower panel; see text for details). The line 
with a slope corresponding to a dissolution timescale of 
2.5 Gyr is displayed in both panels. 

We have also studied the dependence of the cumula- 
tive age distribution on the cluster spatial positions. The 
lower panel of Fig. shows the case for the clusters lo- 
cated at R gc >10 kpc (filled circles - 32 objects) and at 
R gc <10 kpc (open squares - 39 objects), respectively. 
The age distribution for clusters with R gc >10 kpc has 
been rescaled to have the same total number of objects 
as for R gc <10 kpc. The two distributions of points are 
very similar (the result is independent of the definition of 
the Rgc ranges). As a guideline, the best fit single slope 
corresponding to a dissolution timescale of 2.5 Gyr is also 
displayed. This similarity results from the lack of correla- 
tion between age distribution and R gc for the full sample; 
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there is no trend of age with respect to i? gc , with a large 
age spread at any value of the galactocentric distance. All 
of this, in turn, suggests that the cluster formation and 
destruction processes are apparently not correlated with 
the galactocentric distance. 

The upper panel of Fig. [7| shows the cumulative age 
distribution for clusters in two selected ranges of height 
above the Galactic plane \z\. Open squares are clusters 
with \z\ equal to or lower than 300 pc (39 objects), whilst 
filled circles represent clusters at higher distances from 
the plane of the Galaxy (32 objects). Again, the two age 
distributions have been rescaled to the same number of 
objects. It is evident that in this case the two cumula- 
tive functions are different, i.e. the clusters closer to the 
plane follow a relationship much closer to the linear slope, 
corresponding to a dissolution timescale of 2.5 Gyr. The 
more distant ones show a clear excess of clusters in the 
range between 2-4 and 6 Gyr. This is different from the 
conclusion we drew from the distribution for the whole 
sample (Fig. [BJ, the inconclusiveness found there might 
come from mixing two different subsamples. Analogous re- 
sults are found when changing the limits to 250 pc or 350 
pc. The difference is significant (at more than 3<r level), 
even considering the error bars we obtain with a Monte 
Carlo simulation similar to the one discussed for the whole 
sample. This result hints at a more homogeneous creation- 
destruction processe for the clusters closer to the Galactic 
plane, than for their more distant counterparts. However, 
one has to take into account the possibility that this dif- 
ference is, at least partially, an artifact due to the possible 
incompleteness of the OC sample at low \z\ (as discussed 
in JP94), and/or to the cluster orbital motions. The anal- 
ysis by Carraro & Chiosi (1994b) of the orbits of 5 old OCs 
seem to indicate that the observed \z\ for the old OCs do 
not reflect their initial values, due to the rapid oscillatory 
motions of the clusters across the disk. 

3.1. Age-[Fe/H] correlation 

The determination of the age-metallicity relationship for 
halo and disk objects has been the subject of numer- 
ous studies, because it poses a constraint to the chemi- 
cal evolution history of the Galaxy (e.g., Twarog 1980, 
Edvardsson et al. 1993, Carraro & Chiosi 1994, Friel 1995 
and references therein). A first glance at Fig. 01 does 
not show any trend of the cluster age with respect to 
[Fe/H]. A more detailed analysis, however, needs to take 
into account the radial abundance gradient present in 
the Galactic disk (e.g. Friel 1995 and references therein). 
Therefore, we have first determined the relationship - if 
any - between cluster metallicity and galactocentric dis- 
tance; we restricted our analysis to the sample of 38 clus- 
ters with [Fe/H] on the homogeneous scale by GOO - dis- 
played in Fig. |H1 - and we fitted to the data a linear re- 
lationship weighting the various points according to the 
individual [Fe/H] error, thus obtaining 

[Fe/H] = (-0.055 ± 0.019) R gc + (0.37 ± 0.20) (2) 
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Fig. 8. Relationship between [Fe/H] and galactocentric 
distance for the 38 clusters with metallicities on the GOO 
scale. The solid line displays the linear best-fit considering 
all 38 clusters; the dashed line shows the best-fit excluding 
Be 29. 

with a statistically significant slope, in very good agree- 
ment with the value — 0.059±0.010 dex kpc -1 estimated 
by Friel et al. (2002) in their sample of 39 clusters. 

The data displayed in Fig.[S]also clearly show that the 
value of the slope might be affected by the most distant 
cluster in the sample, i.e. Be 29. If we exclude this cluster 
the slope becomes steeper 

[Fe/H] = (-0.097 ± 0.023) R gc + (0.77 ± 0.23) (3) 

but still within the range of independent determinations 
(see, e.g., the discussion by Friel et al. 2002). 

We also tried to assess a possible dependence of this 
slope on the cluster ages, as found by Friel et al. (2002). 
As a first test, we divided the sample into two age 
ranges, i.e. clusters with t < 3 Gyr and clusters older 
than 3 Gyr. For the first group we found a gradient 
A[Fe/H]/Ai? gc = -0.13 ±0.02 dex kpc" 1 , and for the sec- 
ond one A[Fe/H]/Ai? gc = -0.08 ±0.03 dex kpc" 1 . These 
two values are different at about 1.5cr, but the higher age 
group shows a flatter gradient. This is the opposite as 
found by Friel et al. (2002) who determined a steeper gra- 
dient for clusters older than 3 Gyr, with respect to younger 
objects. This could be due to the different [Fe/H] scale 
they used, although their different cluster ages may also 
play a role. 

The precise value of the gradient for the older group 
is again affected by Be 29; if we neglect this cluster the 
gradient becomes — 0.024±0.052 dex kpc" 1 . 

This analysis clearly underscores the need to increase 
future sampling of clusters with [Fe/H] on an homoge- 



10 



M. Salaris et al.: Age of the oldest Open Clusters 




-0.50 



0.00 



0.50 



10.0 



5.0 




-2.00 



[Fe/H] B 



-1.00 
[Fe/H] 



0.00 



Fig. 9. Age distribution as a function of [Fe/H] corr ([Fe/H] Fig. 10. Age distribution as a function of [Fe/H] for the 
corrected to the solar circle, see text for details) for the clusters analyzed in this study (filled circles) and the GCs 



38 clusters with [Fe/H] on the GOO scale. 



studied by SW02 (open squares). 



neous scale, in order to conclusively determine the depen- 
dence of the [Fe/H] radial gradient on age. On the other 
hand, models of Galactic chemical evolution do not pro- 
vide a conclusive prediction about the age dependence of 
the radial [Fe/H] gradient. As shown by Tosi (1996), var- 
ious authors predict a [Fe/H] radial gradient which can 
stay constant with time, increase or decrease (e.g. Fig. 5 
in Tosi 1996). 

We also checked the possible existence of an [Fe/H] 
gradient with respect to the height to the Galactic plane; 
we considered the [Fe/H] values the clusters would dis- 
play at the solar galactocentric distance (assumed to be 
equal to 8.5 kpc), by applying the radial gradient correc- 
tion A[Fe/H]/Ai? gc =-0.055 dex kpc" 1 . No statistical sig- 
nificant trend is found. This result is the same as found by 
Carraro & Chiosi (1994a) in a smaller sample of old OCs; 
according to Carraro & Chiosi (1994b), this may be ex- 
plained in terms of rapid oscillatory motions of the clusters 
across the Galactic plane, which tend to erase any preex- 
isting gradient. On the other hand, the same authors con- 
clude that the observed radial abundance gradients should 
not be seriously affected by the orbital motions, at least in 
the limit of the small sample of cluster orbits (5 clusters) 
they analyzed. 

Based on these results, we have again considered the 
individual [Fe/H] values of our 38 OCs with GOO esti- 
mates, corrected to the solar circle ([Fc/H] corr ). The rela- 
tionship between cluster age and [Fe/H] corr is displayed in 
Fig. |21 and it does not show any statistically significant 
trend between these two quantities, just a large age spread 
at any metallicity. This is consistent with the OC results 
by JP94, Carraro & Chiosi (1994a), Friel et al. (2002), 



and also with the findings by Edvardsson et al. (1993) in 
a sample of field disk stars (although other studies find a 
significant age-metallicity relationship in field disk stars, 
e.g. Twarog 1980). 

4. Comparison between GC and old OC ages 

A comparison between the ages of the oldest OCs and the 
GCs in the thick disk and halo provides vital clues to the 
scenario for Galaxy formation. For this comparison to be 
meaningful it is however necessary to ensure that the OC 
and GC ages are on a consistent scale. 

In SW02 we have accurately and homogeneously deter- 
mined the ages of a large sample of 55 GCs, using stellar 
models computed with the same code and the same input 
physics as the models used in this paper. The differences 
with respect to this work are the age dating method and 
the different He enrichment law assumed in the model 
computation, as discussed previously. This latter point 
has been addressed by redetermining the age of the SW02 
sample using models computed with the same He enrich- 
ment law used for the OCs. The net effect is to cause an 
average decrease by 0.7 Gyr of the GC ages with respect 
to SW02 results. 

As for the different age dating methods applied to the 
two samples, we have one cluster - 47 Tuc - in common, 
whose age has been derived with both techniques. A com- 
parison of its age given in Table 1 with the values from 
SW02 corrected for the new He abundances, provides a 
difference of only 0.5 Gyr (SW02 age being lower), well 
within the error bars associated to the individual deter- 
minations. 
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Fig. 11. Histogram of the ages of the GC sample studied 
by SW02, compared to the OC ages derived in this study 
(shaded histogram). 



Figure ^| shows the distribution of ages as a function 
of the observed [Fe/H] for the SW02 sample (with metal- 
licities according to the Carretta & Gratton 1997 scale) 
and the old OCs of Table 1, the corresponding number 
counting as a function of age are displayed in Figure 1111 
The cut-off in the age distribution of GCs at ~12 Gyr 
and the overlap between the tails of the distribution of 
GC and OC ages is evident. The youngest GCs, supposed 
to have been accreted by our Galaxy, have ages compa- 
rable to the ages of the old OCs. This means that these 
accretion processes were acting well after the formation 
of the Galactic disk. NGC 6791 and Be 17, the two old- 
est OCs, have formally the same age as thick disk GCs 
like 47 Tuc and M 71, implying an approximately coeval 
formation for both thin and thick disk. By comparing the 
ages of the oldest OCs with the oldest GCs one derives a 
difference between the start of the formation of the halo 
and of the thin disk of ~2 Gyr . 



5. Summary and conclusions 

In this paper we have extended our previous age determi- 
nations for GCs to old OCs belonging to the Galactic thin 
disk, using as before a morphological age indicator. In case 
of the old OCs, it is the so-called 8(V) parameter defined 
by JP94. We derived a new and homogeneous calibration 
of the <5(V)-£-[Fe/H] relationship from a subsample of 10 
clusters with accurate and deep photometry, [Fe/H] and 
reddening estimates. Distances to these calibrating clus- 
ters have been determined by means of the MS fitting 
technique using field stars with Hipparcos parallaxes, and 



the ages obtained from fitting the appropriate isochrone 
to the absolute brightness of the cluster turn-off region. 

To obtain reliable distances and age determinations 
for the calibrating clusters, a necessary prerequisite is the 
use of consistent metallicity scales for both field stars 
and calibrating OCs. The metallicities for the unevolved 
Hipparcos field dwarfs used in the OC age calibration were 
derived by PSK02, and shown to be consistent with the 
GOO metallicity scale for OCs. Comparison between the 
[Fe/H] estimates by GOO and the recent work by Friel et 
al. (2002) revealed systematic differences, the most ex- 
treme case being NGC 6791, for which the GOO estimate is 
[Fe/H] = 0.40±0.06, whereas Friel et al. found 0.11 ±0.10. 
A further independent check for the internal consistency 
of our distances (and ages) is possible, by requiring that 
the distance to this cluster derived from the MS fitting, is 
the same when using the (B — V) or the (V — I) colour. 
The metallicity dependence of the MS colour is different 
for these two indices (see, e.g. PSK02), therefore the con- 
sistency of the distances obtained with (B — V) and (V—I) 
is a good test for the adopted metallicity scales. 

Keeping the cluster [Fe/H] as a free parameter, we 
found that consistent distances are obtained - irrespective 
of the choice of the cluster reddening - only when [Fe/H] 
is equal to 0.4, or at least not lower than 0.3, i.e. when it 
is homogeneous with the scale used for the dwarfs. With 
our field dwarf [Fe/H] scale a metallicity, e.g., [Fe/H]=0.2 
for NGC 6791 would cause a discrepancy by 0.11 mag be- 
tween the distances inferred from the (B — V) and (V — I) 
colours. 

With the 5(V)-t-[Fe/H] relation given in Eq. (1), we 
then derived age estimates for a total of 71 OCs. Their 
age scale can be merged with the one we published previ- 
ously for 55 GCs (SW02), 47 Tuc (whose age obtained in 
this paper agrees with the one estimated using SW02 tech- 
nique) being the link connecting the two samples. Due to 
our method, the use consistent isochrones and an homoge- 
neous metallicity scale, we not only obtained the first large 
and homogeneous sample of OC ages, but even a reliable 
age scale on which both cluster types can be placed. This 
allows the investigation of questions related to the forma- 
tion of the various components of the Galaxy, halo, thick 
and thin disk. The bulge still awaits investigation, mainly 
due to the problem of strong and differential reddening of 
the bulge cluster CMDs. 

Using the whole GC and old OC sample fFig. 110(1 . we 
determine a delay by 2. Oil. 5 Gyr between the start of the 
halo and thin disk formation. We estimated this value by 
determining the average age (with error) of the two oldest 
OCs (NGC 6791 and Be 17) - formally coeval - which has 
been then compared with the age of the oldest metal poor 
GCs. Liu & Chaboyer (2000) have estimated 2.8 ±1.6 Gyr 
for this time delay, whereas Carraro et al. (1999) found the 
thin disk younger than the halo by about 2-3 Gyr, an age 
difference shorter than the 3-5 Gyr gap determined by 
Sandage et al. (2003). 

We also find that thin and thick disk started to form 
approximately at the same time, since the age of the thick 
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disk globulars is the same, within errors, as the age of 
NGC 6791 and Be 17. 

The age of the oldest OCs is of the order of 10 Gyr, 
compatible with that of the oldest thin disk white dwarfs 
as estimated from the white dwarf luminosity func- 
tion of the solar neighbourhood, which is, according to 
Hansen (1999), between 6 and 11 Gyr. This rather large 
age range depends on the uncertainties in the observa- 
tional data and white dwarf (surface and core) chemical 
compositions; there are also additional uncertainties due 
to the white dwarf model physics (e.g. Salaris et al. 2000). 

Figure |U| clearly demonstrates the absence of any 
age-metallicity relation, consistent with earlier results by 
Carraro & Chiosi (1994a) and JP94. The overall slope of 
the relationship between [Fe/H] and R gc is consistent with 
recent determinations by, e.g., Friel et al. (2002); however, 
we find a decrease of this slope for increasing cluster ages, 
which is just the opposite of the results by Friel et al. 
(2002). We do not detect any correlation between [Fe/H] 
and height above the Galactic plane \z\, nor between age 
and R gc (as in Carraro & Chiosi 1994 and JP94). 

The cumulative age distribution for the full OC sample 
shows a departure from the predictions of constant for- 
mation rate and exponentially declining dissolution rate 
(with timescale of 2.5 Gyr) at the 2a level. No correlation 
between the cumulative age distribution and R gc is found; 
however, there is a significant excess of clusters in the age 
range between 2-4 and 6 Gyr for the population located 
at hig \z\ values, with respect to their counterpart closer 
to the Galactic plane. It is not clear if this difference is 
intrinsic - i.e. related to the position of the cluster at its 
birth - or partly an artifact due to incompleteness of the 
sample (which, according to JP94, should preferentially 
affect clusters with lower \z\) and/or to the cluster orbital 
motion. 
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